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Abstract 


A nonperturbative analytic solution of Green ’s function for heavy 
ion transport in a single medium is extended to multilayer transport 
media. This extension is implemented in the Green ’s junction code 
(GRNTRN), which is being validated with laboratory experiments. 
Good agreement is achieved between an iron radiation beam experi- 
ment and GRNTRN only when the interactions of the iron radiation 
beam with several important layers of material located upstream of the 
target within the beamline are considered. 


Introduction 

Astronauts and occupants of high-altitude aircraft are exposed to heavy-ion cosmic radiation that 
originates from the Sun and galactic sources. The shielding and exposure of these individuals are 
controlled by the transport properties of the radiation through the specific materials. Efficient space- 
radiation-transport computer codes (refs. 1, 2, and 3) have been developed and applied to a wide range 
of missions, but the results of these codes could not be validated accurately in laboratory experiments 
(ref. 4). Recently, the use of a nonperturbative approach to Green’s function solutions has led to an effi- 
cient laboratory code, namely, Green’s function transport (GRNTRN) code (ref. 5), which will be fur- 
ther developed for space-radiation-transport calculations and validated with laboratory experiments. 

Preliminary validation of the laboratory Green’s function code (GRNTRN) has been made by com- 
paring the theoretical results with recent iron radiation beam measurements made at Lawrence Berkeley 
Laboratory Bevalac accelerator by Shinn et al. (ref. 6). The iron beam had been accelerated to 
600 MeV/amu and passed through a series of beam-transport elements, triggering devices, and a lead- 
foil beam spreader prior to emerging from the beam tube and striking a layer of aluminum or polyethyl- 
ene targets. Overall, fairly good agreement between measurements and predictions from GRNTRN was 
reached when the nuclear database used in GRNTRN was improved in the physical formulation and 
corrected for mass and charge conservation (refs. 7 and 8). Furthermore, as the code was modified to 
include the interactions of beam and lead foil, better agreement was reached except that the predicted 
fragmentation was still low by 30 percent for the aluminum target. The code modification was made to 
extend the nonperturbative solutions for a single medium to inhomogeneous, two-layer transport media 
(for lead foil and target material). In this paper, we present further improvement to our analytical results 
by modifying GRNTRN to include interactions from other significant layers of material from the beam 
tube. The development of the nonperturbative method to inhomogeneous, multilayer (three or more lay- 
ers) transport media will be given with a brief review of the analytic properties of the perturbative solu- 
tion of the heavy-ion-propagation problem and a nonperturbative approximation to the full perturbation 
series. 


Transport Equation 

We restrict our attention to the multiple charged ions for which the Boltzmann equation may be 
reduced (ref. 9) to 



3E Sj(E) + Gj 


= yc jk ^ k (x,E) 


k 


( 1 ) 


where §j(x,E) is the flux of ion type j at x with energy E (in MeV/amu), Sj(E) is the change in E per unit 
distance, <3j is the total macroscopic reaction cross section, and <5^ is the macroscopic cross section for 



collision of ion type k to produce an ion of type j. The solution to equation (1) is to be found subject to 
the boundary condition 


Q>j(0,E) = fj(E) (2) 

The corresponding Green’s function is introduced as a solution of 

I - ^ SjW + E, £') = ^ jk G k Jx, E, E') (3) 

k 

subject to the boundary condition that G ■ (0, E, E') = & -b(E - E'). 

jrn Jin 

The preceding equations can be simplified by transforming the energy into the residual range as 


and defining new field functions as 


'> = i, 


E dE' 


0 ~Sj(E') 


(4) 


= 5/F)^.(x,F) 

(5) 


(6) 

fj(rj) = ~Sj(E)^(0,E) 

(7) 


Equation (3) becomes 




with the boundary condition 


( 8 ) 


r j< O = 5 jm d ( r j 


O 


and with the solution to the ion fields given by 


(9) 


V/(*> r j) = E J” G jm( x ’ r f r ' m )fm(r'J dr' m 

m 

2 

Note that V . is the range scale factor such that V -r . = v r and can be taken as v ■ = Z ./A .. 

J J J rn 171 J J J 

Perturbative Green’s Function 

The solution to equation (8) is written as a perturbation series: 

G jm (*» r ? r m ) = r f r m ) 


( 10 ) 


(ID 


where 



and the first collision term is 


G jL > ( x ’ r J ’ r m) = -rzT ex p <-vrvJ 

m j 


where G\ l J(x, r r' ) takes on nonzero values for 

JtJl J 171 


v, 

V 
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V . 

J -(r j + x)<r’ m <-^r j + x 
v m 


(13) 


(14) 


and 


x j = ( V m x + V /j - V m r J / ( V m ~ V j) ( 15 ) 

x m = ( V m r m ~ V “ V /j^ V m ~ V j ) ( 16 > 

Note that in equations (13)— (16), jc • + x m = x for all r . and r\ The significance of x m is that it is the 
distance ion m traveled from the boundary to the collision site at which the ion j was produced and must 
now travel distance Xj before reaching x. Similar results hold for higher order terms (refs. 10 and 11). 

Whenever r.,r' < x, it is relatively easy to show that 

J • *1 


!r r J ’ r ’"‘ )dr j = V a m) (17) 

where the integral limits are the interval limits in equation (14) and g(j, m ) is the g-function of two argu- 
ments defined by the following relations: 

g(j) = exp (-Ojx) (18) 


and 


gU V j 2 ’ ‘"Jn’jn+0 ~ 


g(j\’j 2 ’ Jn-l’Jn) jn-l*Jn+0 


°jn + 1 in 


(19) 


The first collision Green’s function term may be approximated by the spectral averaged value as 


G^. l \x r . r )~ i r r -'“^(l)/ w r 

^jm p r m> r _ r J KJ jm r p r m> ar 

ju jl Jr jl 


_ V j°jm8U ’ m ) 


xV m~ V j 


( 20 ) 


and is functionally dependent only on depth of penetration x. Its contribution applies only over the 
energy range of equation (14) and is null outside. The higher order terms are approximated by their 
spectral averages as 


G jm( X ’ O’ r 'm ) 


I 
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( 21 ) 
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Note that Gj2 (*, r j, r' m ) is purely dependent on x for i > 0, and we represent the expression as Gj2 (■*) 
(ref. 4). In terms of equations (1 1)— (21), the solution to equation (1) becomes (ref. 4) 

¥/(*» r j) = e xp {-OjX)fj(rj + x) + X G j2( x ) 

m, i 

In equation (22), r mu and r m j are given by the upper and lower limits of the r , satisfying 
equation (14). The symbol F m (r m ) refers to the integral spectrum 

Fm(r' m ) = f , fm(r) dr (23) 


Fm(r’,)-F m (r‘) 


( 22 ) 


Nonperturbative Green’s Function 

We now introduce nonperturbation terms for the summation in equation (22). First, we note that 


Sj m M = & Jm g( m '> + a jmgU,m) + Y d ^j k ^ km g(j,k,m)+... 

k 


(24) 


A term-by-term comparison of the series (22) with the series (9) shows that, in the spectral average 
approximation, we may write the approximate solution as 

V v /[«j m (*)-exp(-Cj:) 6 , ] r . , -I 

V, (*>'•;) = ex P (-V j x)f(r j + x) + Y J - J - _ T | V _ v i [ F ’n(r mu )-F m (r ml )^ (25) 

m \ m J\ 

which is a relatively simple quantity once gj m (x) is evaluated (ref. 4). The advantage of equation (23) 
is that gj m (x) satisfies the convolution product relation (ref. 12) as 


gjmW = 'Lg i k(*~y)gkJy) (26) 

k 

for any positive values of x and y. A solution for small x is easily obtained from the first few terms of 
equation (24). Equation (26) may then be used to propagate the function g -(x) over the entire solution 

jin 

space, taking arbitrarily large steps. 

Nonperturbative Method for a Shielded Medium 

The major simplification in the Green’s function method results from the fact that the scaled spec- 
tral distribution of secondary ions to a first approximation depends only on the depth of penetration, as 
seen in equations (20), (21), and (22). Our first approach to a multilayered Green’s function will rely on 
this observation and assume its validity for multilayered shields. 

Consider a domain labeled as 1 that is shielded by a second domain labeled as 2; the number of type 
j ions at depth x in 1 due to type m ions incident on domain 2 of thickness y is 

«12 jm^y) = X SU*WS 2 * m (>’) (27) 

k 

The leading term in equation (24) is the penetrating primaries as 

= ex P(-°l j x - Vljy YOjm + lgl2jm( x > - v ) “ ex P(-d ly j: - <J 2j y)O jm \ (28) 
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where all higher order terms are in the bracket of equation (28). The first term of the scaled Green’s 
function is then 


G \2jm( x <y> r j> r 'm') = e xP(-G\jX-G 2 jy'> & jm & l x+r J-( r 'm-Py)] < 29 ) 

where p is the range scale factor for the two media 

p = Ry{E)IR 2j {E) (30) 

which is fairly constant above 10 MeV/amu or so. We take a single value for p that corresponds to 
600 MeV/amu. The secondary contribution is similarly found from equation (14) as 

v ■ v . 

- l (r j + x + py)<r' m <-±r j + x+py (31) 

v m v m 

from which the average spectrum is evaluated. The full approximate Green’s function is then 


G \2jm( x ’ >> r j’ r ’m> “ ex P(-°l/ - a 2/) 8 > 8 ( Jt + M + r j ~ r ’m> 

, y) - e x p(-qj/ - 

( x + py)|v ra -v y 


(32) 


Equation (32) is our first approximation of Green’s function in a shielded medium of two layers, and is 
easily modified to multiple layers. (See appendix.) 

For the first spectral modification, the first collision term has the properties 


V ) 0 Um eXp( “°lm A: “ a 2m>’) 


G 12jm( x ’y’ r j’ r 'm ) = •! 


V - V ■ 
m j 
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(33) 


These properties are used to correct the average spectrum to 


G 1 2jm^y’ r r r 'n) = 


V WjnMyy +b ( )(r . 

(x + py)\v m ~Vj\ jm{ ,y { m m 


0) 


where g]2j m ( x > y) * s the first collision term of equation (32) and 


(34) 


>* -L. y 

mu ml 


(35) 


is the midpoint of r' m between its limits given by equation (31). The bj m term of equation (34) has the 
property that 



which ensures that the first term of equation (34) is indeed the average spectrum as required. The spec- 
tral slope parameter is 


, , _ v 7 V m[°ljm ex P(- a lm' c - 0 2m3’)-Cr 2j „,exp(-a 1 ^-0 2 ^)] 

im(X ’ y) ~ (^py)(v m -v.)|v m -v.| (37) 

A similarly simple spectral correction could be made to the higher order terms. The spectral correction 
given in equation (37) is included in the present Green’s function code. 

Spectra for Laboratory Beams 

The boundary condition appropriate for laboratory beams is given by 


/..(£) = ——E exp 
J J2k g 


~-(E-E 0 ) 

2a 2 




(38) 


where Eq is the mean energy of beam with an energy spread a. The cumulative spectrum is given by 


If ( E - £ 0 yi 


J2<5 yj 


(39) 


The cumulative energy moment needed to evaluate the spectral correction is 


‘ r( E ~ E 0T 

a 

■ (E-Eq) 2 - 


+ exp 

,J2k 

. 2a 2 _ 


(40) 


The average energy on any subinterval (Ep E 2 ) is then 

£;(£,) -g/g 2 ) 

F/£,)-F.(£ 2 ) 

The beam generated flux is 


(41) 


¥/(*» ^ r j) = exp(- G lj x - G 2j y)}j( r j + x + py) 

+ 'L G jm y)lFm(r' mu ) ~ Fm(r' m ,)] 
m, i 

+ y) l^(M) ~ r'J[F m (r mu ) - F m {r’ ml )] (42) 

m 

where E is evaluated with equation (41), with E ^ and E 2 as the lower and upper limits, respectively, 

associated with r', and r' . 

ml mu 

The linear energy transfer (LET) distributions behind aluminum and polyethylene targets were 
measured with a stack of four CR-39 plastic nuclear track detectors (PNTD’s) placed behind the target 
(ref. 13). The incident beam intensity was monitored by a thin layer of CR-39 PNTD placed immedi- 
ately in front of the target. Calibration of these detectors with various ion beams of known LET resulted 
in a response function that is approximately Gaussian with a LET-dependent width. The detectors and 
targets were run in good geometry so that acceptance corrections are not required. The detector response 
function is used to compare the experiment with the calculated LET distributions from GRNTRN, 
assuming 0.2 percent beam-energy spread in equation (38). 
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r\ 

Figures 1(a) and (b) show measured and calculated LET spectra behind 2-g/cm aluminum and 

r\ 

4.6-g/cm polyethylene targets, respectively. The calculations include the interactions between the iron 
beam and 2.24 g/cm of the “thin” lead-scattering foil placed upstream of the target. The energy of the 
iron beam prior to its interaction with the two-layer lead-target combination was 557 MeV/amu. The 
calculated LET spectra for fragments are in good agreement for the polyethylene target experiment 
(fig. 1(b)), but are 30 percent lower than the measurement for the aluminum target (fig. 1(a)). Note that 
the spread in each isotope peak results mostly from the detector-response function rather than the beam- 
energy spread. 

Some of the materials used in the beamline, other than lead foil, that could cause significant modifi- 
cations to the beam prior to its interaction with the target were further identified. In our latest calcula- 
tion, a 0.16-cm-thick (0.43-g/cm 2 ) aluminum exit plate and two polystyrene scintillators of 0.16-cm and 
0.32-cm thickness (total 0.49 g/cm 2 ) were included in our transport analyses. Equations (42) and (32) 
were modified for a four-layer (lead foil, aluminum exit plate, polystyrene, and target) configuration 
and implemented in GRNTRN. The improved result for the 2-g/cm aluminum target appears in 
figure 2. For the polyethylene target, the results of a four-layer calculation (not shown) were not appre- 
ciably different from those of the two-layer calculation, as was expected. Because the fragmentation for 
the polyethylene target is far greater than those for the aluminum target, the contributions from the alu- 
minum exit plate and polystyrene are far less significant than those from the lead-foil and polyethylene 
target. The beam energy for the four-layer shield-target combination was inferred, taking into account 
the energy loss through these layers, to be 583 MeV/amu. The latter energy value is much closer to the 
beam energy at extraction (600 MeV/amu), which indicates the four-layer configuration is a good 
approximation to the exact one. 

Concluding Remarks 

The iron radiation beam experiment with CR-39 plastic nuclear track detectors (PNTD’s) provides a 
fairly accurate laboratory validation of Green’s function transport (GRNTRN) code. Through the multi- 
layer analyses, it is concluded that for thinner targets, the effect of lead-scattering foil, beam-transport 
elements, and triggering devices in modifying the characteristics of the beam is rather significant. 
Laboratory validation up to this point provides confidence in the validity of the Green’s function 
method and, in an approximate sense, in the correctness of the nuclear physics theory. Any refinement 
in the nuclear database used in GRNTRN ultimately will depend on better resolution of each isotope 
measurement with a state-of-the-art detector system (Zeitlin et al., Radiat. Meas., vol. 23, no. 1, 1994) 
and more sophisticated analyses, including multilayer configuration. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
December 16, 1996 
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Appendix 

Derivation of Approximate Green’s Function for Three-Layer Configuration 

The formalism presented in the preceding text is extended to a three-layer configuration as follows. 
The solution to equation (26) in a three-layered medium is 


Smjm( x ’y< z '> = 'Lsi jk Mg2ki(y)s 3 , m (z) 
kl 

The leading term is the penetrating primaries, and equation (Al) may be written as 
Snijm( x <y<z.) = exp (-<jj jX - o 2 jy - a 3 jz)8 jm 

+ ig 123 jJ x ’ y- Z > ~ ex P(°l j x ~ °2 jy ~ °3 /) 5 jml 
The scaled Green’s function is then 

G 123 jm( x ’ y> z > r p r 'm) * ~ ex P (“°l/ “ ^ ~ a 3j z > S jm 8 ( x + P2 ? + P3* + r j ~ r 'm ) 

t v /gl23 jm( x ’ * g) - exp(-g ly .x - c 2j y - g 3 ;z)5 y - J 
(x + p 2 y + p 3 z)( v m - v J .) 


(Al) 


(A2) 


(A3) 


where p 2 = R 1 j(E)/R 2 j(E) and p 3 = R^j(E)/R 2 -(F). The range condition of equation (17) 
becomes 


rAO + x + p 2 y + P 3 *) - r m S Z Lr j + * + P 2 y + P 3 * 

v m v m 

The spectral corrections are similarly derived. 


(A4) 
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(a) 2.24 g/cm 2 Pb + 2 g/cm 2 Al. 



Figure 1. Transport measurements and calculations for lead and target (two-layer) combination for 
557-MeV/amu iron beam. 



Figure 2. Transport measurements and calculations for four-layer shield and target combination 
(2.24 g/cm 2 Pb -I- 0.42 g/cm 2 Al + 0.49 g/cm 2 polystyrene + 2 g/cm 2 Al) for 583-MeV/amu iron 
beam. 
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